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Abstract. We study the generation and maintenance of large-scale magnetic fields in barred galaxies. We take a 
velocity field (with strong noncircular components) from a published gas dynamical simulation of Athanassoula 
(1992), and use this as input to a galactic dynamo calculation. Our work is largely motivated by recent high quality 
VLA radio observations of the barred galaxy NGC 1097, and we compare our results in detail with the regular 
magnetic fields deduced from these observations. We are able to reproduce most of the conspicuous large-scale 
features of the observed regular field, including the field structure in the central regions, by using a simple mean- 
field dynamo model in which the intensity of interstellar turbulence (more precisely, the turbulent diffusivity) is 
enhanced by a factor of 2-6 in the dust lanes and near the circumnuclear ring. We argue that magnetic fields can 
be dynamically important, and therefore should be included in models of gas flow in barred galaxies. 
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1. Introduction 

Barred galaxies possess remarkably rich and unusual mag- 
netic structures. The first systematic observations of po- 
larized radio emission from a sample of barred galaxies 
(Beck et al. 2001a, hereafter Paper I; see also Beck et al. 
1999) have revealed widespread polarized emission in the 
bar region, indicative of magnetic fields ordered at scales 
of order 1 kpc. The best examples to date are the nearby 
prototypical barred galaxies NGC 1097 (Beck et al. 1999, 
see Fig. ||) and NGC 1365 (Paper I; see also Fig. 7 in 
Lindblad 1999), and these galaxies show broadly similar 
magnetic structures. 

The field in the bar region, obtained from the B- 
vector of the polarized radio emission at centimeter wave- 
lengths where Faraday rotation is small, appears approx- 
imately to be aligned with the gas streamlines result- 
ing from gas dynamical simulations (Athanassoula 1992; 
henceforth A92); notably, it is the velocity field in the 
reference frame corotating with the bar that exhibits the 
alignment. Such an alignment is perhaps not surprising in 
view of the strong gas velocity shear (Beck et al. 1999). In 
NGC 1365 (Lindblad et al. 1996) and NGC 1530 (Reynaud 
& Downes 1998) the gas has strong velocity gradients near 



Send offprint requests to: D. Moss 



the dust lanes, consistent with the predictions of gas dy- 
namical models. For NGC 1097, velocity data with suf- 
ficient resolution are not available. The total nonthermal 
radio emission, a tracer of the total (regular-|-turbulent) 
magnetic field, exhibits ridges near the dust lanes (see also 
Ondrechen & van der Hulst 1983) that occur because of 
the strong gas compression downstream of the bar's ma- 
jor axis (A92). Thus, the hydrodynamic shock appears to 
have a counterpart in the total magnetic field (but not in 
the regular field - see below). 

There are, however, a number of puzzling features in 
the magnetic field which are especially evident in the regu- 
lar magnetic field (as opposed to the turbulent field whose 
scale is of order 100 pc). Firstly, the contrast in the to- 
tal magnetic field strength near the dust lanes is signifi- 
cantly lower (about 2 in NGC 1097, see Beck et al. 1999) 
than that in gas density (about 20 or more), as inferred 
from models (A92, Englmaier & Gerhard 1997) and from 
CO observations (e.g. for NGC 1530, Reynaud & Downes 
1998). This suggests that the total magnetic field strength 
B does not scale with the local gas density p, as would fol- 
low from a one-dimensional shock compression in the dust 
lanes, B (x p, and even B oc p^" is too strong a depen- 
dence. Beck et al. (1999) argue that the low magnetic field 
compression could be due to enhanced turbulent diffusion 
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Fig. 1. The radio continuum emission in the southern part 
of NGC 1097, with the galactic centre at top right, ob- 
served at A6.2 cm with the VLA (see Paper III for details). 
Contours give the total intensity (representing the total 
magnetic field), dashes show the observed polarization 
plane rotated by 90° (representing the regular magnetic 
field uncorrected for Faraday rotation). The Faraday rota- 
tion is about ±20° at this wavelength (although it can be 
larger in localized regions). The length of the dashes is pro- 
portional to the polarized intensity. The half-power width 
of the synthesized beam (shown in the lower right cor- 
ner) is 10". The background optical B-band image (grey 
scale) was obtained at the Cerro Tololo Interamerican 
Observatory by H. Arp. The bar major axis is at an an- 
gle —30° and the disc major axis (the line of nodes) is at 
—45°, both with respect to the south- north direction. The 
south-western side of the galaxy is closer to us. 



downstream from the shock, but this idea needs quantita- 
tive verification. 

Moreover, the polarization intensity (a measure of the 
strength of the regular magnetic field) in NGC 1097 has a 
broad local maximum upstream of the dust lanes, where 
gas density has deep minimum. The regular magnetic field 
shows hardly any enhancement in the dust lanes. It is also 
surprising that the regular magnetic field does not show 
any discontinuity associated with the shock. The sharp 
turn in magnetic lines leading to a depolarized strip de- 
tected by Beck et al. (1999) has been resolved in later 
observations (Beck et al. 2001b, Paper III) into a smooth 



turning of the magnetic lines in the region 1 kpc upstream 
of the dust lanes (Fig. ^). The turn is even smoother in 
NGC 1365 and occurs about 3 kpc upstream of the com- 
pression region (Paper III). Thus, there is no clear signa- 
ture of any shock-related compression or sharp deflection 
in the regular magnetic field. It appears that shock effects 
in the magnetic field may have been mostly eradicated by 
the enhanced turbulent diffusion in the sheared shock re- 
gion. The verification and quantification of this effect is 
one of the main goals of this paper. 

There are thin, feather-like, dust filaments bent in a 
similar way upstream of the main dust lane, seen in the 
optical photographs of NGC 1097 and NGC 1365. The ob- 
served magnetic field is well aligned with these filaments, 
especially at a distance of 1-3 kpc upstream of the shock 
front in NGC 1097 (see Fig. |l]). Similar filaments have 
been observed in NGC 5383 by Sheth et al. (2000) who 
observe that they are associated with H ii regions and pro- 
pose that the high density and low shear in the filaments 
make them hosts to star forming regions. If the dust fil- 
aments trace the streamlines, the velocity vectors must 
deviate from the streamlines predicted for this region by 
gas dynamical simulations (see e.g. Fig. 12 in Roberts et 
al. 1979, Fig. 2b in A92, and Figs 16 and 26 in Lindblad 
et al. 1996 where the streamlines have sharp cusp rather 
than turn smoothly upstream of the shock front). If the 
magnetic field is sufficiently strong in that region it may 
control the shape of the dust filaments; however, the oc- 
currence of a dynamically dominant magnetic field has to 
be confirmed observationally and then explained. 

These facts imply that the physics involved is more 
complicated than that of a simple shock with passive mag- 
netic field. The equipartition strength (with respect to 
the cosmic rays) of the observed regular magnetic field 
is about 7 /iG upstream of the dust lane in the southern 
part of NGC 1097 (Beck et al. 1999). Taking a gas num- 
ber density of n = 0.1-0.2 cm~'^ (from A92), we obtain an 
Alfven speed of 30-50 km s~^. This is much larger than the 
speed of sound (and, presumably, the turbulent velocity if 
it has the standard value of lOkms^^) and is comparable 
to the shearing velocity in most of the bar region. Hence 
the magnetic energy density seems to be comparable to 
the kinetic energy density in the regular shearing motion 
and much exceeds the thermal energy density. It follows 
that a fully consistent model of gas flow in barred galaxies 
should include magnetic flelds, a factor systematically ne- 
glected in all gas dynamical models of barred galaxies. In 
this paper we discuss models of regular magnetic fleld evo- 
lution in the velocity fleld resulting from Athanassoula's 
(1992) gas dynamical model in order to assess the extent 
of the problem, i.e. to establish which features of the ob- 
served magnetic field can be readily explained by existing 
knowledge of the gas flow in barred galaxies and which 
appear to demand more advanced MHD modelling. 

Models of magnetic fleld evolution in barred galaxies, 
without the inclusion of dynamo action, have been con- 
sidered by Otmianowska-Mazur et al. (1997), who used 
a velocity fleld resulting from A^-body simulations (with 
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38,000 stars and 19,000 molecular clouds). As usual with 
A^-body simulations with a modest number of particles, 
their density distribution has only a moderate enhance- 
ment near the bar major axis where dust lanes occur in 
real galaxies. The simulated magnetic field concentrates 
in the spiral arms outside the bar region and is weak in 
the bar. Since the magnetic field is not supported by any 
dynamo action, it decays in less than 5 x 10^ yr. 

Dynamo models for barred galaxies have been consid- 
ered by Moss et al. (1998, 1999) who used the gas velocity 
field obtained from A^-body simulations with 200,000 stars 
and 40,000 gas particles. Their density field shows appro- 
priately stronger density contrasts and the inclusion of 
dynamo action results in persistent magnetic structures. 
However, the distribution of particles in such models is 
sparse and uneven, so it is difficult to define reliably a 
velocity field in low-density regions. This leads to serious 
complications in magnetic field modelling, and so we pre- 
fer to use a velocity field obtained from gas dynamical 
simulations. 

Our model focuses on the regular magnetic field, and 
the effects of interstellar turbulence on the magnetic field 
are parameterized in terms of turbulent transport coefh- 
cients, including a standard a-effect (Moffatt 1978, Parker 
1979, Krause & Radler 1980). The formulation would be 
similar for alternative mechanisms, such as a magnetic 
buoyancy driven dynamo, as suggested by Moss et al. 
(1999b). The key feature of our model is the inclusion 
of a realistic representation at high spatial resolution of 
the velocity field in barred galaxies. 



We also report briefly in Sect. 4.5 on our discovery of 



a rather remarkable dynamo effect in galactic discs with 
positive dynamo number, where the dynamo is driven by 
noncircular motions rather than by differential rotation. 
This may be of interest in the wider context of dynamo 
theory. 

2. The alignment of magnetic and velocity fields 

As noted by Beck et al. (1999), the regular magnetic field 
in the bar region of NGC 1097 seems to be aligned with the 
velocity field of A92, especially in regions with stronger ve- 
locity shear. This can be seen from a comparison of Fig. m, 
where we present the polarization map of NGC 1097, with 
the magnetic field orientation indicated with dashes, and 
Fig. g where a velocity field from A92 is shown. Such an 
alignment is not typical of normal spiral galaxies where 
magnetic field lines are inclined to the streamlines by 10°- 
30°, presumably due to the dynamo action (e.g.. Beck et 
al. 1996, Beck 2000, Shukurov 2000). In the presence of a 
strongly sheared velocity, the local structure of the mag- 
netic field will be controlled by the local velocity shear. 
In barred galaxies, the shear of the noncircular velocity 
field is strong enough to make the form of their magnetic 
fields markedly different from those found in normal spiral 
galaxies. 

Ignoring the effects of dynamo action, the large-scale 
field will be frozen into the flow in regions with i?,n — 



uL/rj ^ 1, where u ~ lOOkms"^ is the regular shearing 
velocity, L ~ 3 kpc is its scale in the bar region, and i?,„ 
is the magnetic Reynolds number based on the turbulent 
magnetic diffusivity "q. Thus, the field will be aligned with 
the flow if 77 ^ 10^^ cm^ s~^. 

However, the alignment of magnetic field and the flow 
can be affected by dynamo action even at large values 
of i?,n. Dynamo action is needed to maintain the global 
magnetic field against the effects of winding by differen- 
tial rotation and tangling by turbulence, which would lead 
eventually to enhanced Ohmic decay. Therefore, we also 
require that the dynamo is unable to misalign the field 
and the streamlines: the local growth rate of the magnetic 
field 7 must be smaller than the shear rate, i.e. 7 ^ u/ L 
with 7 ~ D'^/'^rj/h'^, where D — auh^ /{Lif) is the local 
dynamo number and a is the alpha-coefficient. This yields 
aL/uh ^ 1, or Lvt/hu ^ 1, where vt is the turbulent ve- 
locity and a < Vt (since a cannot exceed the turbulent 
speed - see, e.g.. Sect. V.4 in Ruzmaikin et al. 1988). In 
normal galaxies where u is dominated by the streaming 
velocity produced by the spiral pattern and so L/h 2± 2-3 
and ut/u ~ 1, this inequality is not satisfied and we can 
not expect strong alignment between the streamlines and 
magnetic lines. Indeed, the magnetic pitch angle (i.e. the 
angle between the regular magnetic field and the stream- 
lines) is about 1/3 radian, plausibly consistent with this 
estimate. On the other hand, the shear rate u/L is sig- 
nificantly larger in barred galaxies and we can expect a 
much closer alignment in the regions where L ^ 3 kpc and 
u Jl 60kms~^, where the above inequality is satisfied. In 
other words, we expect that barred galaxies contain large 
regions where dynamo action is overwhelmed by the local 
velocity shear resulting in a tight alignment of the mag- 
netic field with the shearing velocity (more precisely, with 
the principal axis of the rate of strain tensor). On the 
other hand, there may be regions of enhanced diffusivity 
and/or reduced shear where the alignment is reduced (see 
Fig. I). 

The velocity field, but not the magnetic field, looks 
different in the inertial and corotating frames. We can ex- 
pect a rather close general alignment between B and u 
in the reference frame corotating with the bar for the fol- 
lowing reason. Nonaxisymmetric magnetic field patterns 
must rotate rigidly to avoid winding up by differential ro- 
tation (and it is the dynamo which can maintain such 
fields). With a nonaxisymmetric perturbation from the 
bar (or spiral arms), the magnetic modes that are coro- 
tating with the perturbation will be preferentially excited 
(e.g. Mestel k Subramanian 1991; Moss 1996, 1998; Rohde 
et al. 1999). Thus, the regular magnetic field will corotate 
(or nearly corotate) with the bar, and this is a physically 
distinguished reference frame. (All magnetic field configu- 
rations discussed in this paper, except those in Sect. 4.5, 
do exactly corotate with the bar.) So, an approximate 
alignment between B and u is expected in the corotat- 
ing, but not in the inertial, frame. 
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Fig. 2. The velocity and gas density (grey scale, lighter 
shades correspond to higher values) fields in the frame 
corotating with the bar, from Model 001 of A92. The sense 
of rotation is clockwise. The bar major axis is at an angle 
of 45° to the vertical, the frame radius is 8 kpc and the 
corotation radius is re « 6 kpc. 

3. The model 

3.1. The gas velocity field 

The velocity field was taken from data supplied by 
E. Athanassoula, corresponding to Model 001 shown in 
Fig. 2 of A92. This model is two-dimensional, with u — 
{ux{x, y), Uy{x, y), 0) in cartesian coordinates {x, y, z). The 
gas velocity and density fields are steady in the frame ro- 
tating with the bar (the corotating frame). The stellar 
bar extends approximately between the ends of the dust 
lanes (shown as a light shade in Fig. H). We reproduce in 
Fig. g the velocity field in this frame rotating in a clock- 
wise sense with angular velocity of about 34 km s~^ kpc~ , 
placing the corotation radius at about r = re ~ 6 kpc. 
The streaming velocities relative to the rotating frame 
are as large as lOOkms"^ near the shock fronts, which 
are slightly offset from the major axis of the bar. Regions 
where the sheared velocity is comparable with the rota- 
tional velocity are widespread in the bar region. A shock 
near the dust lanes is just resolved in the data. The region 
in which the velocities are available is 16 kpc across. 

The magnetic fields obtained using this velocity field 
will be compared with observations of NGC 1097 where 
the corotation radius is Vc ~ 12 kpc (assuming a distance 
to NGC 1097 of 17 Mpc); the length scales are there- 
fore renormalized correspondingly before any comparison 
is made. However, when discussing the computational re- 
sults in Sects. Q and 15 we retain the original length scale 
of A92 with corotation at 6 kpc. Since the model of A92 



refers to a generic barred galaxy rather than to the galaxy 
NGC 1097 specifically, our comparison with observations 
can only be qualitative. 

3.2. The dynamo model 

The basic dynamo model, as applied to galaxies with 
strong noncircular motions, is described in Moss et al. 
(1999); it uses the 'no-z' approximation (Moss 1995), re- 
placing derivatives perpendicular to the galactic midplane 
(z-direction) with inverse powers of /i, the disc semi- 
thickness. This is in some ways consistent with the two- 
dimensional flow model of A92. The mean field dynamo 
equations are solved for the magnetic field components in 
the directions orthogonal to z. [B^ is given in principle 
by the condition V ■ B = 0.) The resulting magnetic field 
can be thought of as an approximation to the field in the 
mid-plane, or to represent values averaged vertically over 
the disc. 

We made one significant amendment to the model. 
Phillips (2001) has showed that the dynamo growth rates 
in the 'no-z' model can be made to agree more closely 
with those of the asymptotic analysis (e.g. Ruzmaikin et 
al. 1988) by including a correction factor of 7r^/4 into the 
terms representing diffusion in the vertical direction, and 
so we also included these factors. Thus the local marginal 
dynamo numbers for our models can now be expected to 
be fairly directly comparable to standard values of about 
10. 

The large velocity shear present in barred galaxies re- 
sulted in unsatisfactory numerical behaviour of the nu- 
merical algorithm previously employed (e.g. Moss et al. 
1998, 1999; Moss 1995), and so a new version in carte- 
sian coordinates was written, using a second order Runge- 
Kutta method for time-stepping, and second-order accu- 
rate space discretisation. This code solves the 'no-z' dy- 
namo equations in the frame corotating with the bar. The 
earlier calculations were carried out using the auj formula- 
tion, but now we allow for the regeneration of both merid- 
ional and azimuthal regular magnetic field by interstellar 
turbulence (the a^Lo dynamo). The dynamo equation has 
the form 



iV77 xB + aB)- Vx (7/VxB),(l) 



dB „ , 
— -=Vx MXB 
ot 



where u is the regular velocity field described in Sect. |3.1| , 
a parameterizes the dynamo action of the interstellar tur- 
bulence, and 77 is the turbulent magnetic diffusivity (see, 
e.g., Ruzmaikin et al. 1988, Beck et al. 1996). In Eq. (|I]), 
we ignore any vertical (z-wise) variation of 77, but allow a 
variation in the disc plane. (In principle we might, in the 
spirit of the no-z approximation, have approximated the 
effects of a halo diffusivity larger than that in the disc by 
an estimate of the form drj/dz ^ ?7(x, y)/h.) The equation 
is solved in a cylindrical region with radius 8 kpc. We al- 
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low for nonlinear dynamo effects leading to the saturation 
of the magnetic field growth by adopting 



ao 



a{r) 



1 + S,B'^/AT:p{r)v, 



(2) 



where ao (a typical amplitude of a) is a constant, a(r) 
is the normalized background value, and p and Vt are the 
gas density and turbulent velocity. We are thus assuming 
that the large-scale magnetic field significantly reduces the 
a-effect when its energy density approaches that of the 
turbulence; ^ is a constant, introduced to suggest some of 
the uncertainty about the details of this feedback. When 
estimating field strength from our models we set ^ = 1, 
but this fundamental uncertainty should be kept in mind. 
(Indeed, it is quite possible that, rather than being a con- 
stant, ^ could depend on local conditions.) We assume that 
Vt = lOkms""'^, and p is taken from the same gas dynami- 
cal model of A92 as the velocity field u. It is often assumed 
that a oc oj with u the angular velocity of rotation. We 
first consider models with a constant apart from quench- 
ing effects, i.e. a{r) = 1, and then put a(r) — uj{r)/uj{r^), 
where we arbitrarily choose r^ = 0.3i?, and R — 8kpc. 
(Note that ui{r^)/Vti^ « 0.22 whe re the unit angular veloc- 
ity r^o is introduced in Sect. 3.3.) We choose this normal- 
ization so as to have the dimensionless quantities of order 
unity in the outer parts of the galaxy. 

We can estimate the z-averaged vertical field in the 
disc, by using the condition V • B = 0, and estimating 
dBz/dz by Bz/h. This vertical magnetic field (not shown 
in Fig. and other similar figures), is small on average, 
being about 10 times weaker than the average horizontal 
field. However, the vertical field is significant in regions 
with pronounced structure in the horizontal field (e.g. the 
shock front and the central region) where it can be com- 
parable to the horizontal field. 

Note that the dynamo equations with the quadratic 
nonlinearity (0) allow the transformation B -^ ~B\ there- 
fore the magnetic field vectors in the figures shown below 
can be reversed. 

Equation (|l|) is nondimensionalized in terms of length 
i?, the characteristic disc radius, and the magnetic diffu- 
sion time across the galactic disc, h? /rjQ with h the disc 
scale height and 770 a typical value of ry. Magnetic field 
is measured in units of Bq — ^J An p^^v^ / ^ (corresponding 
to equipartition between the turbulent and magnetic en- 
ergies), where po is the maximum density. The induction 
effects arising from turbulence and large-scale motions are 
quantified by dimensionless dynamo numbers 



Ra 






Ru. 



noh^ 



where subscript zero denotes a characteristic value of the 
corresponding variable. From here onwards we refer only 
to dimensionless quantities, unless explicitly stated other- 
wise. 

Although we are studying an a^uj dynamo character- 
ized by two separate dimensionless parameters Ra and 
i?[^, it is still useful to consider the dynamo number 



D = —RaRu as a crude measure of the dynamo intensity. 
i?a) is positive here because it is defined in terms of a typ- 
ical angular velocity flo > 0. In the standard asymptotic 
analysis for galactic dynamos (see Sect. VI. 4 in Ruzmaikin 
et al. 1988), D is defined in terms of rduj/dr, which is 
negative. Thus we have introduced a minus sign into our 
definition of D so that, consistent with Ruzmaikin et al., 
the standard galactic dynamo has D < 0. Our results are 
primarily a function of D (for given /„, /^, see Eq. (ph) 
being quite insensitive to the relative values of Ra and i?cj 
in the parameter range considered; this indicates that ro- 
tational velocity shear dominates in the production of the 
azimuthal magnetic field (i.e. we have approximately an 
auj dynamo) . The dynamo can maintain the regular mag- 
netic field if l^^l > £'cri where D^r is a critical dynamo 
number, which depends weakly on details of the model 
(Sect. VI. 4 in Ruzmaikin et al. 1988); as a rough esti 
mate, Dcr — 8 in the simple model discussed in Sect. 4.1 



This marginal value (obtained with the velocity field of 
A92) is close to that found for a calculation using only 
the axisymmetric part of the azimuthal velocities. 

The computational domain formally is —1 < x,y < 
+1, but we solve only in the region r = {x^ + y'^Y^'^ < 1, 
with boundary conditions B^, By = Q a.t r = I. Field 
amplitudes are found to be small near the boundaries of 
the domain in most cases considered, and so the bound- 
ary conditions remain self-consistent throughout the sim- 
ulations. The region r > 1 is ignored; formally B — 
there. Our choice is conservative, being arguably the least 
favourable for dynamo action; for example dynamo action 
occurring outside r = 8 kpc could feed field into the region, 
interior to this radius, where the equations are solved. 

Our standard computational grid has 160 mesh lines 
in both the x- and j/-directions, uniformly distributed over 
~1 ^ x,y < 1; we made a trial integration on a finer grid, 
finding insignificant differences in our results. 

3.3. Units 

We choose the size of the formal computational domain to 
be 16 kpc, consistent with that of A92. Thus R — 8 kpc, 
and we choose h = 0.4 kpc, although the overall nature of 
the solutions is relatively insensitive to the latter choice: h 
appears in the dimensionless control parameters Ra and 
Rui , and the overall field strength depends on their product 
D. With r^o = 150 km s^^ kpc~^ (of order the value in the 
inner Ikpc radius), we have Rui « 72/7726, where 7726 = 
770/10^^ cm^ s~^ and, if ao is a few kms~-^, then Ra will be 
of order unity. We assume that Vt = lOkms"^ everywhere, 
so that, with po « 3.5 x 10~^^gcm~'^, the unit magnetic 
field strength is Bq = 6.6 x 10-^^"^^^ G. 

Neither vt nor Cg are well known from observations of 
barred galaxies. However, Englmaier & Gerhard (1997) 
found that sound speeds less than about 20kms~^ are 
necessary to produce shocks shifted away from the bar 
axis. We note in this connection that random velocities of 
molecular clouds in the central regions of NGC 1097 are as 
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Fig. 3. Contours of 77/770 — 1 wit h /^ 
u from A92, as used in Sect. 4.2 



3 in Eq. {}^ with 
The 5 contour levels are 



equally spaced with a range from 0.10 to 0.90. As in Fig. g. 
the frame radius is 1.7 times the corotation radius. 



high as 35kms~^ (Gerin et al. 1988), i.e. 3-4 times larger 
than the turbulent velocities at larger distances from the 
centre. A similar central enhancement in the velocity dis- 
persion of molecular gas has been detected in NGC 3504 
(Kenney et al. 1993). We take Vt = lOkms"^ as a uniform 
background value, but the turbulent intensity is assumed 
to be enhanced in the shock fronts and in the central re- 



gion as described in Sect. 3.4 



A representative value of the turbulent magnetic dif- 
fusivity is 770 — hvtl ~ 10^^cm^s~^ for v^ = lOkms"^, 
where / ~ 100 pc is the turbulent scale. A convenient esti- 
mate of a is (Sect. V.4 in Ruzmaikin et al. 1988) 



a ~ mm | — — , Vt 
h 



0.5kms 



n 



20kms-ikpc" 



lOkms" 



(3) 



,(4) 



where we neglect possible local enhancements of turbulent 
velocity beyond lOkms""'^. With the above values of I, h 
and ilo, we obtain ao — 4kms~^. Thus, our representative 
values of the dynamo coefficients are Ra = 5 and Rui = 
70. Since Eq. (g) is just an order of magnitude estimate, 
we consider that a variation of a factor 3 in the dynamo 
coefficients is quite acceptable. 



3.4. Enhancement of turbulence by shear 

We allow for the possibility that both 77 and a may be 
enhanced where shear flow instabilities are likely to occur, 
by writing 

a = ao (l-l-/„-^^ ), 77 = 77o(r) (n-/^- ), (5) 

\ '-'max / \ '-'max / 

where /„ and /^ are constants and S — \dux/dy\ + 
\duy/dx\, u = (uxjUy) and S'max is the maximum value of 
S. For given values of Ra, R^, and/or D, the local dynamo 
action is weakened if /r; > 0, and enhanced if /q, > 0. 
Contours of 77/770 — 1 with /^ = 3 are shown in Fig. Bt 
our models have magnetic diffusivity enhanced in the dust 
lanes and in the innermost 1-2 kpc. 

The effects of varying a and 77 cannot be disentangled 
completely since it is the ratio a/77^ that affects magnetic 
fields generated by an au) dynamo. Therefore, the result of 
enhanced a can be reproduced, quite closely, by reducing 
77 instead. Thus, in order to keep the models as simple as 
possible, we preferred to take /q. = in most simulations. 
Furthermore, with S defined after Eq. (||), the turbulent 
transport coefficients a and 77 with fa, frj 7^ would be 
enhanced even in rigidly rotating regions. The most im- 
portant result of this is that 77 can be unrealistically large 
in the central parts of the disc in models discussed below. 
A perhaps more physically meaningful prescription could 
he S = \dux/dy + ti>| + \duy/dx — lo\. However, the effect 
of the latter refinement is reduced when working in the 
rotating frame and, moreover, it would plausibly be simi- 
lar to that of reducing /^. So we leave this refinement for 
future, more advanced models. A trial calculation showed 
that the effect is small when /,, ;$ 3 and ao is constant. 

We suggest that the enhancement of the turbulent dif- 
fusivity in the shear flow (the dust lane) is due to the de- 
velopment of flow instabilities. The enhancement will be 
sufficient to reduce the contrast in the magnetic field if the 
turbulent diffusion time is shorter than the advection time 
across the dust lane width: d? /r] < d/cg with d ~ 0.5 kpc 
in NGC 1097 (Regan et al. 1995, 1997) and we have as- 
sumed that the transverse velocity immediately behind 
the shock is equal to the speed of sound Cg — lOkms"-^ 
(e.g. Roberts et al. 1979, Englmaier & Gerhard 1997). This 
yields 77 > 1.5 x lO^'^cm^s"^ for the turbulent magnetic 
diffusivity within the dust lanes, so 77 would have to be 
increased about tenfold from the background value 770, in 
order to reduce the field contrast. 

The Kelvin-Helmholtz instability is an obvious can- 
didate for turbulence amplification in the dust lanes (see 
Townsend (1976 ) or Terry (2000) for a recent discussion 
of turbulence in shear flows and Ryu et al. (2000) and 
Briiggen & Hillebrandt (2001) for a discussion of mixing 
enhancement by the instability). The growth rate of the 
instability in a piecewise continuous flovif is given by 

1/2 



^™ = Yd"" 



2K- 



1-2K- 



, sinh K 



expK 

where K = kd with k the wavenumber in the direction 
of the shear flow, d ~ 0.5 kpc is the transverse width 
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of the sheared region and U ~ 100 kms"^ is the ve- 
locity difference across it. The growth rate has a max- 
imum, 7max ~ O.lU/d^ for K w 0.8. The resulting en- 
hancement in the turbulent velocity can be estimated as 
Vt ~ i'toexp(7i„ax'i/cs), where wto ~ Cs ~ lOkms"^ is the 
turbulent velocity upstream of the dust lane, and we have 
assumed that the transverse velocity immediately after the 
shock is equal to the speed of sound, so the residence time 
of the matter in the unstable region is given by d/cg. This 
yields 7maxrf/cs — 1, so the turbulent velocity in the dust 
lane is expected to be 2-3 times larger than the upstream 
value of lOkms^^. The most unstable mode has a scale 
of I = 2TTd/0.8 w 4kpc. The resulting enhancement in the 
turbulent diffusivity rj ci ^vtl can be as large as by a factor 
of 100, but this enhancement will be limited by the fact 
that only a part of the flow produced by the instability will 
be randomized, and simultaneously the effective value of 
I will be reduced. Arguably, a better estimate of 77 result- 
ing from the instability is 77 ~ ji't 7max — Swtoc' — ISwto^j 
implying enhancement by a factor of 50. The above esti- 
mates are upper limits, referring to the maximum rather 
than the mean enhancement. The enhancement of 77 used 
in our models ranges up to a factor of /^ -I- 1 = 6, and the 
case fjj — 3 is illustrated in Fig. 0. 

A magnetic field parallel to the shear velocity can 
suppress the Kelvin-Helmholtz instability if Va > U 
(Chandrasekhar 1981), but this inequality is not satisfied 
in our case. 

Apart from the Kelvin Helmholtz instablity, 
enhanced small-scale three-dimensional motions 
can arise from local perturbations to the gravi- 
tational field as discussed by Otmianovifska-Mazur 
et al. (2001). 



4. Results 

4.1. Our basic model 

We first examined the simplest dynamo model with 
a(r) = 1, fa = fri = 0, i.e. 77 is constant, unmodified by 
the shear, and a is a function only of the magnetic field 
and gas density, via the a-quenching. We experimented 
with several parameter combinations, and first discuss 
here a model with 770 = 2 x 10^^ cm-^ s~^ so R^ = 36. With 
Ra = 0.3 (giving D = —10.8), the dynamo is slightly su- 
percritical in the central parts of the galaxy (r < 1-2 kpc) 
and subcritical at larger radii. The eventual magnetic field 
configuration is steady in the corotating frame, and its 
structure is shown in Fig. 0, superimposed on the gas den- 
sity. Here, and in the other figures, we have smoothed the 
magnetic field over a scale of about 0.6 kpc, to approxi- 
mate the resolution of the observations. For comparison, 
the radio observations have a beamwidth of 10" (0.8 kpc) 
for NGC1097 and 15" (2 kpc) for NGC 1365. We see that 
the magnetic field is strongly concentrated to the central 
regions (r ^ 2 kpc in this model with uniform diffusiv- 
ity, refiecting the strong dynamo action in this region 
where the angular velocity and its shear are maximal. 



This feature remains even when the dynamo number is 
very substantially increased. In this model, the maximum 
field strength of about 0.5i?o is reached at r « 0.7 kpc 
(corresponding to about 3/iG); the field is negligible be- 
yond r ~ 2 kpc. Observations of NGC 1097 suggest that 
the ratio of the regular field strength in the inner region 
{r/vc — 0.1) to that in the outer bar region is about 2, 
and so this model is unsatisfactory. Note that this feature 
is little altered by the smoothing to the resolution of the 
observations. 

Even when \D\ is increased by an order of magnitude, 
the model magnetic field remains strongly concentrated 
to the central regions. Note that in all our models the 
field strength increases with \D\, but in order to obtain 
field strengths comparable with those observed near to 
the corotation radius, it is necessary to increase \D\ to 
values at the margin of what is plausible, and the field 
in the inner regions is then significantly larger than ob- 
served. The overall field structure is relatively insensitive 
to the magnitude of D. An obvious problem with this ba- 
sic model with i] = const is that it yields an unrealistically 
weak magnetic field in the outer regions and perhaps a too 
strong contrast in B between the dust lane regions and 
those upstream of them. Thus we now relax the condition 
that 7/ be constant. 

4.2. Models with turbulent diffusivity depending on 
velocity shear: uniform a 

An obvious refinement of the model is to consider the tur- 
bulent magnetic diffusivity to be modulated by the shear 
rate, i.e. to put /^ > in Eq. (^), keeping the background 
value uniform, i.e. rio{r) = const. Note that, as f,j and 
thus the magnetic field dissipation increase, so does the 
marginal dynamo number, for example Dcr ~ —27 when 
/,, = 5, and even D = — 11 is just subcritical when /,, = 1. 
Also, broadly speaking the effects of taking /q > for 
given Ra and of increasing Ra with fa — are similar, 
and so we keep /„ — 0. We now set rjo ~ 10^® cm^s"^, 
so that Ri^ — 72 and adopt Ra = 3. Thus the models we 
discuss now are all substantially supercritical, with the lo- 
cal values of D (calculated from the local values of 51) of 
order —(10^-10'^) at r ^ 3 kpc and about —15 near the 
corotation radius. 

In Fig. u we show the field structures for calculations 
with fri = 1, 3, 5. Again, the magnetic fields are steady in 
the rotating frame. Increasing /^ clearly reduces the dom- 
inance of the central field. Also, increasing 770 to 2 x 10^^ 
cm^s"^ {Rui = 36) produces little overall effects on the 
field structure, but does affects the field strength. When, 
for example, R^ = 36, fr, = I, Ra = 0.6 {D = —22.5), 
the maximum field strength is about 0.7Bo] this figure in- 
creases to about 6.5i?o when R^^ — 72, /^ = 1, Ra = 3 
(D = —216, Fig. ^) Of this increase, a factor of about 
4 can be attributed to enhanced induction by rotational 



shear which yields B ~ Bq {\D\/Dc 
to the noncircular velocities. 



1) , and the rest 
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Fig. 5. Vectors of the horizontal magnetic field for the 
models with 5(r) = 1, /„ = 0, Ra = 3, i?cj = 72. (a) 
/,, — 1, (b) fri = 3, and (c) /,, = 5. The circumscribing 
circle has radius 8kpc. Shades of grey show gas density, 
with lighter shades corresponding to higher values. The 
magnetic field vectors have been rescaled as S — > B/B"" 
with n = 0.3, to improve the visibility of the field vectors 
in regions with small B. 



A feature of all these solutions is that there is a broad 
field minimum in the bar region, but with strong mag- 
netic ridges at the positions of the dust lanes and field 
enhancement in the central part where the circumnuclear 
ring is observed. This agrees well with the overall dis- 
tribution of polarized intensity in NGC 1097 and 1365. 
Nevertheless, the local magnetic energy density near the 
major axis (upstream of the dust lanes) exceeds the tur- 
bulent kinetic energy density. As we discuss in Sect. 5.1 



both the model and observed field strengths have a well 
pronounced structure within the bar region in addition to 
the ridges elongated with the dust lanes and the central 
ring. 



There is a relatively strong, quasi-azimuthal field up- 
stream of the shock, especially near the ends of the bar. 
This feature also agrees well with the observations, cf. 
Fig. 1^. For the dynamo parameters of the models illus- 
trated in Fig. pi the dynamo is locally supercritical nearly 
everywhere. However, the field structure is rather insensi- 
tive to the value of D. To illustrate this, we show in Fig. |^ 
the field vectors with only a sHghtly supercritical dynamo 
number [rjo ^ 2 x lO^^cm^s'S i?„ = 0.6, R^; = 36). The 
field configuration is remarkably similar to that of Fig. pk, 
with a much larger dynamo number, except perhaps near 
the outer boundary. This supports the idea that the field 
in regions upstream of the dust lanes at larger galacto- 
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centric distances is mainly produced not by local dynamo 
action but rather by advection over almost 180° from re- 
gions at small radii (where the dynamo action is stronger) 
on the other side of the bar. We note that the half-rotation 
time of the gas at r ~ 5 kpc is shorter than magnetic dif- 
fusion time h?/riQ ~ 5 x 10^ yr. If the dynamo is strong 
enough {Ra ^ 3 for R^ — T2,) the advection produces a 
magnetic configuration where the magnetic energy den- 
sity exceeds that of the interstellar turbulence in broad 
regions, including the region upstream of the shock front. 
The steep rotation curve at small radii and the large local 
dynamo numbers found there suggest that the dynamo 
is particularly efhcient in the inner region. This field is 
then transported outwards by the noncircular velocities 
and further enhanced by shear. 

Immediately downstream of the shock position, the ve- 
locity and magnetic fields are closely aligned. The align- 
ment is fairly good throughout the galactic disc, with the 
mean angle x between magnetic and velocity vectors being 
about X « 20° ((cosx) ~ 0.93) in the model of Fig. |a. 
As can be seen in Fig. R, the alignment is especially close 
outside the bar and slightly reduced within the bar and 
at the corotation radius. The reduction in the bar can be 
plausibly attributed to enhanced magnetic diffusivity in 
the dust lanes and that at corotation to the small local 
values of u in the rotating frame. 



4.3. Models with turbulent diffusivity depending on 
velocity shear: 5 ex a; 

We now discuss models in which we set a oc uj{r) (see 
Sect. 3.2), whilst keeping /„ ~ 0. Since the dynamo action 
now becomes more intense in the disc centre, the field 
strength is more strongly concentrated to the centre for 
given /^, than when a is uniform. As /^ is varied, these 
results can be broadly summarized by saying that the field 
structure is very similar to that when a{r) is uniform, if 
a larger value of /^ is now taken. For example, we show in 
Fig. H the field structure when R^ = 72, /^ = 3, i?Q = 3, 
which can be compared with that shown in Figs. Ha and 
^3, with air) — 1 and /^ = 1,3 respectively. We note that 
the effect of setting a (x Lu(r) is the opposite of increasing 
fri - the former enhances the central field concentration, 
the latter decreases it. But the field structure outside of 
the inner 1-2 kpc is altered only slightly by such changes. 



For this model, the maximum field strength is about 
3.9i3o, attained at r/vc ~ 0.1 (see also Fig. O). We show in 
Fig. S a grey scale plot of B^ for the model of Fig. @ where 
magnetic field enhancements in the dust lanes, the spiral 
arms and the central region are prominent. The field is es- 
pecially strong in the circumnuclear region and at the ends 
of the bar. Trailing magnetic arms, emanating approxi- 
mately from the ends of the bar, are a common feature 
of the modelled and observed magnetic structures. (This 
feature is also seen in the model for the weakly barred 
galaxy IC 4214 of Moss et al. 1999.) 



This class of models exhibits the best agreement with 
observations and we discuss them in detail in Sect. 0. 

4.4. Magnetic fields in the absence of dynamo action: 
the case of pure advection 

In this section we attempt to elucidate the roles of the two 
key ingredients of our model, stretching by the regular ve- 
locity and the a-effect. In Fig. |l^ we show the field struc- 
ture for a computation with R^ — 72, /,, = 5, and R^ = 
(no dynamo action) . In the absence of dynamo action this 
field decays exponentially, with an e-folding time of about 
10^ yr. This decay time scale is about 5 times shorter than 
in normal galaxies because of the enhanced velocity shear 
that reduces the field scale and so facilitates its decay. 
Since most barred galaxies in the sample of Paper I do 
exhibit detectable polarized radio emission and thus pos- 
sess regular magnetic fields, we conclude that the need for 
dynamo action in barred galaxies is even stronger than in 
normal galaxies. 

Looking at the spatial structure of the decaying mag- 
netic field, we see that some of the features of the dynamo 
maintained fields shown in Figs. |^ and ^ are present, par- 
ticularly the field minimum upstream of the shock, and the 
trailing magnetic arms outside of the corotation radius. 
However there is no hint of the relatively strong fields at 
small radii seen in all models with dynamo action present. 
In addition, the relative magnetic field strength near the 
bar major axis is now too low. 

Thus it appears that the field structure in the outer 
regions is relatively insensitive to the parameters of the 
dynamo model, but that dynamo action is essential to 
maintain the fields over time scales of gigayears and to 
produce strong magnetic fields in the region of the cir- 
cumnuclear ring and upstream of the shock front. 

4.5. Results with positive dynamo number 

We made some experiments with Ra < 0, i.e. positive 
dynamo number D. This is motivated by suggestions that 
this may be the case in accretion discs (Brandenburg et al. 
1995) even though the origin of turbulence there is quite 
different from that in galaxies. It is a familiar result for 
axisymmetric galactic discs that the marginal positive dy- 
namo numbers for dynamo excitation are much larger in 
magnitude than the marginal values for the usually con- 
sidered case with D < - i.e. dynamos with D > are 
much harder to excite than those with D < 0. As men- 



tioned in Sect. 3.2, in the models with D < the marginal 



values of D for the models with the full A92 velocity field 
are essentially unchanged from those found by including 
only the axisymmetric rotational velocities. When D > 0, 
the situation is quite different. A dynamo is excited by the 
full velocity field at a marginal dynamo number of mag- 
nitude considerably less than required to excite the corre- 
sponding dynamo driven by the circular motions alone. In 
other words, when D > the dynamo is driven by shear 
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due to noncircular velocities rather than from differential 
rotation. 

We illustrate our results by discussing a simulation 
with R^ — 72, fff = 5, fa = 0, a{r) = 1 and Ra — —6, a 
somewhat supercritical value. The magnetic fields are now 
oscillatory, with a dimensionless period of about 0.48, cor- 
responding to about 1.1 X 10^ yr. We see in Fig. |ll] that 
the regions of strong magnetic field correspond to the re- 
gions of strong velocity shear. The overall field structure 
varies significantly during the cycle period, and at certain 
times the amount of structure in the magnetic field is much 
greater than at others. The magnetic field strength has a 
maximum in the inner bar region and in two elongated fea- 
tures parallel to the shock fronts. These structures persist 
through nearly all of the oscillation period. However, the 
field geometry within the structures changes with time as 
a magnetic field reversal along a line approximately par- 
allel to the shock front propagates clockwise in Fig. |ll|. 
In the lower left-hand part of Fig. |ll|, weak field directed 
approximately radially outwards is visible. At a slightly 
later time, this has grown in strength and this region eats 
in to the adjacent area of approximately inwardly directed 
field, strengthening the field reversal. A little later still, the 
reversal has disappeared, and then the cycle repeats. We 
emphasize that, for this calculation the field decays when 
the nonaxisymmctric velocities are removed. 

It is clear that these dynamos with positive dynamo 
numbers generate magnetic fields that are rather different 
from those observed. 



5. Discussion and comparison with observations 

5.1. The magnetic field in the shocl< region 

We show in Fig. O the variation of B^ along two cuts 
through the shock front region for the model illustrated 
in Fig. pi illustrating the variation in the magnetic field 
strength between the post- and pre-shock regions. Also 
shown is the observed profile of polarized synchrotron in- 
tensity along similarly chosen cuts in NGC 1097. In the 
case of constant energy density of cosmic-ray electrons, the 
polarized intensity depends on B^ , where B± is the mag- 
netic field component perpendicular to the line of sight. 

Both the model magnetic field and the observed polar- 
ized intensity reveal several peaks in B^ on each side of the 
shock front and a pronounced minimum at the shock front 
itself. The model and observations differ in detail, which is 
no surprise as inter alia the velocity field and gas distribu- 
tion used in our models was not designed specifically to fit 
NGC 1097. We find that the best correspondence between 
the B^ variations from the models and the observed po- 
larized intensity profiles is obtained in models with 5 oc w 
(Sect. 4.3), with Ra ^ 1.5 when R^^ = 72. Larger values of 
\D\ give weaker peaks near the bar, but taking too small 
a value of |Z)| gives rather low field strengths. 

Due to the inclination of NGC 1097 (45°), the dif- 
ference between B^ and B^ is typically a factor of 2. 
Projection effects reduce the polarized intensity upstream 



of the shock from where magnetic field has significant com- 
ponent parallel to the line of sight, but not in the dust lane 
where the value of B± is closer to that of B. This will sup- 
press the amplitude of the secondary peak in polarized 
intensity in the pre-shock region. With allowance for the 
unaccounted projection effect, our model shows remark- 
ably good agreement with observations in the number and 
width of the magnetic peaks in the shock region. 

Our models reproduce a smooth turn of the 
magnetic field vector upstream of the dust lanes 
observed in NGC 1097 (and also NGC 1365). This 
turn may be difficult to see in the figures presented 
here because we have shown magnetic vectors only 
at a fraction of mesh points. It can be seen from 
Fig. 1^ that the alignment between the magnetic 
and velocity fields is reduced 12 kpc upstream of 
the shock fronts (dust lanes). Several effects appar- 
ently contribute to smear and displace the turn in 
the magnetic field with respect to that in the ve- 
locity field. We argue in Sect. 4.2 that magnetic 
field can be advected to the pre-shock region from 
smaller radii (the elongated shape of the stream- 
lines responsible for that can be clearly seen in 
Fig. 2b of A92, but not that easily in Fig. ^here). 
The admixture of magnetic field generated else- 
where must smear any sharp structures produced 
locally. Furthermore, magnetic field will tend to be 
aligned with the prinicpal axis of the rate of strain 
tensor rather than the velocity field itself, and so 
its turn can start long before it becomes visible in 
the velocity vectors. 

Other effects can be envisaged that might, in 
principle, result in a smooth turn of magnetic field 
upstream of the shock fronts. For example, this 
can be a manifestation of the dynamical influence 
of the regular magnetic field on the velocity field. 
However, this is hardly the case in our models 
where the regular magnetic field does not reach 
equipartition with even the turbulent kinetic en- 
ergy in the regions concerned (see Fig. g). Another 
possibility is that the turn in the observed field is 
smoothed by the presence of a magnetic halo with 
a regular magnetic field smoother than in the disc; 
however, there is so far no evidence for a polar- 
ized halo in NGC 1097. A further option is that 
the magnetic field is anchored in a warm compo- 
nent of the interstellar medium that has different 
kinematics than the cold gas (i.e., higher speed of 
sound); our preliminary analysis of this possibility 
has given no evidence for this, but we shall return 
to it elsewhere. 



5.2. Radial profiles of magnetic field 

We show in Fig. |l^ the dependence of the azimuthally 
averaged value of B^ on radius for the model shown in 
Fig. y and a similar model with larger /^, and com- 
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Fig. 12. Profiles of the square of the 
horizontal magnetic field along cuts 
perpendicular to the shock front of 
the model shown in Fig. H at (a) 
1.5 kpc and (b) 3kpc galactocentric 
distance. The total length of each 
cut is 10 kpc, approximately sym- 
metrical about the position of the 
shock front, with distance increasing 
in the upstream direction. The cor- 
responding cuts (with scaled lengths 
of 20 kpc; 20" corresponds to r/rc ~ 
0.14) through the map of A6 cm po- 
larized intensity (represented by the 
vector length of Fig. |l|) are shown 
in (c) and (d) for galactocentric dis- 
tances of 3 and 6 kpc. The cuts are 
shown for the same values of r/r^ in 
the model and NGC 1097 and ex- 
tend over the same linear length. 



pare it with the observed dependence of polarized inten- 
sity on radius. The model reproduces the peak observed 
close to the centre and has a realistic radial scale. The 
slower decrease of polarized intensity with galactocentric 
radius in NGC 1097 can be explained by a steeper rota- 
tion curve in the A92 model, which has the turnover ra- 
dius at r/vc « 0.25 whereas it occurs at about r/vc = 0.4 
in NGC 1097 (Ondrechen & van der Hulst 1989). The 
model magnetic field has a pronounced subsidiary peak 
at r/rc ~ 0.7; this peak is related to the strong magnetic 
arms visible in Fig. |9|. A similar peak in the observed radial 
profile is less pronounced and occurs at a smaller fractional 
radius. 

In our models, the regular magnetic field van- 
ishes on the rotation axis because of its symme- 
try properties. However, the observed polarized 
intensity remains significant at r = 0. We discuss 
implications of this in the next section. 

5.3. The azimuthal structure of the global magnetic 
field 

The azimuthal structure of the regular magnetic field is 
a useful diagnostic of both the overall structure of the 
gas velocity and the field generation mechanisms. We 
have fitted the polarization angles observed in NGC 1097 



Table 1. The global magnetic structure in NGC 1097 
represented in terms of the amplitudes of the azimuthal 
Fourier modes TZm, their pitch angles Pm and phases 
(3rn with m the wavenumber. The bottom line shows az- 
imuthally averaged magnetic field strengths obtained as- 
suming He — 0.1 cm^'^ and h — 0.4 kpc. 







Radial range 


(kpc) 






2.9-3.7 


3.7-4.5 


4.5-5.3 


7^o 


radm"^ 


-152 ±6 


-141 ±13 


-152 ±17 


Po 


degrees 


18 ±1 


27 ±2 


22 ±2 


7^l 


radm~^ 


-75 ±8 


-88 ±12 


-90 ±10 


pi 


degrees 


-50 ±13 


-71 ±13 


-58 ±7 


I3i 


degrees 


-119 ±6 


-98 ±6 


-122 ±8 


7^2 


radm~^ 


110!?a5 


57!L5 


75 +r 


P2 


degrees 


Ollao 


_7+o 

' -102 


-31 IL 


132 


degrees 


13 ±4 


11 ±15 


28 ±8 


(B) 


AtG 


5.3 1";-^ 


4.9 ±0.6 


5.3 tf" 



at AA3.5 and 6.2 cm, for three radial rings between 2.9 
and 5.3 kpc with a model where the magnetic field is ex- 
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pandcd into azimuthal Fourier modes. The model and 
the fitting procedure are discussed by Berkhuijsen et al. 
(1997) and Fletcher et al. (2001); Table |I] gives details 
of the fits. The polarization angles have been obtained 
from Stokes parameters averaged over 20° sectors in the 
0.8 kpc wide rings. The results are presented in terms of 
TZm = 0.81 Bmrich, p„i and j3m, where Bm (in /iG) is 
the amplitude of the magnetic mode with azimuthal wave 
number m, ric is the thermal electron density in cm^'', h 
is the scale height of the thermal ionized layer in pc, Pm 
is the pitch angle of the mode m (i.e. the angle between 
the field direction and the circumference; p > indicates 
a trailing spiral) and (3^, is the azimuthal angle, measured 
counter-clockwise from the northern end of the major axis 
of the galaxy, where the mode m (m ^ 0) has maximum 
strength. TZm is related to the Faraday rotation measure; 
TZm = 100 rad m~^ corresponds to a field strength of about 
3 /iG for Tic — 0.1 cm~'^ and h = 400 pc. Note that the lim- 
ited azimuthal resolution of the fits (to < 2) does not allow 
the fitting of the sharp deflection of the B vectors in the 
shock wave region (Fig. |l|). 

The amplitudes of the three lowest azimuthal Fourier 
modes in the observed regular field of NGC 1097 have 
ratios 1:1: 0.5. The corresponding values for the models 
shown in Figs. pk,b,c and Fig. are all 1 : : 0.5. The 
axisymmetric mode is dominant in both the observed field 
and in all these models. The relative amplitude of the 
m = 2 mode for these models is also in agreement with 
the observations, having about half the amplitude of the 
m ~ mode. 

The magnetic fields in our models do not contain the 
m = 1 magnetic mode present in the observed field. This 
is due to the strict symmetry of the model velocity field 
(it does not contain any modes with odd to) and the in- 
ability of the dynamo to maintain the to = 1 mode on 
its own. The only mode actively maintained by the dy- 
namo in our models is the axisymmetric one, to = 0, and 
the m = 2 mode arises via distortion of the basic field 
by a fiow with a strong m — 2 component. In turn, other 
magnetic modes with m even are also maintained. In real 
galaxies such a high degree of symmetry of the velocities 
cannot be expected, and thus a wider range of azimuthal 
modes will be present in the magnetic field. We note that 
NGC 1097 is perturbed by a companion, so the real ve- 
locity field will certainly be more complicated than that 
used in the dynamo simulations. Any m = 1 component 
of velocity will then generate a corresponding component 
of magnetic field from the to = field by the u x B 
term in the induction equation (e.g. Moss 1995), although 
whether the amplitude of the resulting to, = 1 mode is 
large enough to be consistent with the observed value can 
only be determined by a detailed calculation. Naive nu- 
merical experimentation suggests that the amplitude of 
the TO, = 1 velocity component would have to be com- 
parable with that of the m = 2 component if a, m = 1 
magnetic field of the strength implied by Table |l| is to be 
produced solely by the u x B interaction. Anisotropy 
in the a-effect, neglected here, can also excite the 



TO, = 1 magnetic mode via dynamo action (Rohde 
& Elstner 1998). Neglecting the to = 1 mode, the im- 
plication is that all the models presented in Figs. H and H 
reproduce the global azimuthal structure sufficiently well. 

The absence of modes with odd m explains why the 
model magnetic field vanishes at the disc centre as all 
modes with even to, have _B = at r = from symmetry 
considerations. Models with a more realistic velocity field 
can have strong magnetic fields with m — 1,3,... at the 
disc centre. The significant observed polarized intensity 
at r = (see Fig. |l^) may be due to magnetic modes with 
odd TO,. 

The estimates of azimuthally averaged magnetic field 
strength resulting from the fits of Table ^ shown in the 
bottom line of the Table, have been obtained using an 
electron density of rie = 0.1 cm~'^ and an ionized disc 
scale height oi h = 0.4 kpc. (The errors given do not in- 
clude any uncertainty in ric and h.) Both parameters are 
poorly known for barred galaxies. The diffuse Hq, flux in 
barred galaxies is only moderately lower than in normal 
spirals (Rozas et al. 2000) and the Hi scale height is com- 
parable to (or somewhat larger than) in the Milky Way 
(Ryder et al. 1995), so electron densities can be expected 
to be similar in different galaxy types. The electron den- 
sity at comparable fractional radii (with respect to coro- 
tation) of 4-5 kpc in the Milky Way are about 0.1 cm^^ 
(Taylor & Cordes 1993). This justifies our crude estimate 
ric = 0.1 cm~^ for r ~ 5 kpc in NGC 1097. A significantly 
lower value of Hq would imply a systematically larger mag- 
netic field strength than estimates resulting from energy 
equipartition with cosmic rays. 

5.4. The field in the central ring 

In the top right of Fig. Hi the field at a galactocentric 
radius of about 1.5-2 kpc can be seen to have a strongly 
nonaxisymmetric structure apparently dominated by the 
TO, = 2 mode with a significant m = mode also present. 
The overall field structure features strong magnetic fields 
at those azimuthal angles where the dust lanes intersect 
the nuclear ring. This can be compared with the magnetic 
fields in the inner regions of our models (seen most clearly 
in Figs. 0, Ha, H where the field in the inner regions is 
emphasized) . The form of these fields is remarkably similar 
to that shown in Fig. n. We emphasize that the fact that 
the observed polarized intensity does not vanish at small 
galactocentric radii (see Fig. |l3|), as it should for modes 
with even to, indicates that the modes with odd m should 
become increasingly dominant at smaller radii. Angular 
momentum transfer due to the magnetic stress produced 
by the regular field can provide a mass infiow rate into the 
central region of about 1 Mq yr~^, which is comparable to 
that required to feed the active nucleus (Beck et al. 1999). 

6. Conclusions 

Our simulations confirm that conventional mean-field dy- 
namo models can reproduce magnetic field structures that 
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are similar to those observed in barred galaxies. Even 
though we have used a velocity field model for a generic 
barred galaxy (Athanassoula 1992), our magnetic config- 
urations show an overall agreement with observations of 
the prototypical barred galaxy NGC 1097. Our models are 
most successful when they include enhanced turbulence in 
regions of strong shear, most notably near the shock fronts 
offset from the bar major axis. The best overall agreement 
with observations seems to be provided by models with 
a (X uj and /^ « 3, although the overall appearance of the 
field in the outer regions is not very sensitive to the choice 
of parameters. 

We argue that regular magnetic fields can be dynami- 
cally important in large regions within the corotation ra- 
dius. The energy density of the magnetic field in our mod- 
els and, plausibly, in real barred galaxies exceeds that of 
the interstellar turbulence and is largely controlled by the 
local shear in the regular (noncircular) velocity. This is in 
a striking contrast to the situation in normal spiral galax- 
ies where regular magnetic fields and turbulence are close 
to energy equipartition. 

Our models can successfully reproduce some salient 
features of the observed field structures, notably the un- 
expectedly strong field upstream of the shock, the gentle 
bending of the field across the shock, and the overall struc- 
ture of the field near the galactic centre. There are well 
developed trailing spiral arms present outside of the coro- 
tation radius, which begin approximately at the ends of 
the bar. 

The intensity of interstellar turbulence must be en- 
hanced in the dust lanes and circumnuclear region in or- 
der to obtain a radial contrast in the regular magnetic field 
similar to that observed in NGC 1097. We argue that this 
enhancement can arise naturally due to shear flow insta- 
bilities. The enhancement of turbulence may have impor- 
tant implications for star formation and gas dynamical 
models in barred galaxies. Enhancements in turbulent ve- 
locity implied by CO observations for the central parts of 
NGC 1097 (Gerin et al. 1988) and NGC 3504 (Kenney et 
al. 1993) are consistent with those used in our models. 

The initially surprising bending of magnetic lines up- 
stream of the shock can be explained by (i) advection of 
magnetic field from other regions, and (ii) by the align- 
ment of B with the principal axis of the rate of strain 
tensor rather than the velocity itself. Note however that 
advection in an approximately uniform (non-sheared) flow 
does not cause alignment of u and B, even in the perfect 
conductivity limit. 

In conclusion, we emphasize three important findings 
of this paper. 

1. Barred galaxies are different magnetically to 'normal' 
spiral galaxies, in that the field strengths are signif- 
icantly higher, and the alignment between magnetic 
and velocity vectors is, on average, closer. 

2. The magnetic field may be dynamically important: 
Alfven speeds can be comparable with noncircular ve- 
locities (although generally a little lower). Such strong 



fields require relatively, but not unrealistically, strong 
dynamo action [Ra ^ 3 for our models with 770 = 10^^ 
cm^s^-'^). 
3. In order to model satisfactorily the magnetic fields it 
is necessary for the turbulence to be locally enhanced 
in the vicinity of the dust lanes and the circumnuclear 
ring, in a manner consistent with other evidence. This 
will affect other aspects of modelling of these galaxies. 
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Fig. 4. Vectors of the horizontal magnetic field for the 
model with R^ = 0.3, R^i = 36, fa = frt — 0. The cir- 
cumscribing circle has radius 8 kpc. Shades of grey show 
gas density, with lighter shades corresponding to higher 
values. The magnetic field vectors have been rescaled by 
B -^ B/B"^ with n = 0.3, to improve the visibility of 
the field vectors in regions with small B. Panel (a) shows 
the field vectors over the entire computational domain, 
whereas in (b) the inner region has been omitted, in order 
to show the structure in the outer regions of the disc. 
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Fig. 6. As in Fig. ga, but for i?„ = 0.6, R^ = 36, /^ = 1, Fig. 8. As in Fig. |5|a, but for a{r) ex uj{r) and i?„ = 1.5. 
giving weak dynamo action. 




Fig. 9. Grey scale map of B^ (unsmootlicd) for the model 
Fig. 7. The alignment of the regular magnetic and veloc- shown in Fig. pi. Lighter shades indicate higher values, and 
ity fields in the corotating frame is illustrated with level the contours shown are for B^ — 1, so the energy 
contours of cosx — \u ■ B\/uB for the model of Fig. ^. density of the regular magnetic field exceeds that 
The levels shown are 0.4 (dotted), 0.7 (dashed) and 0.95 of turbulence inside the contours. The circumscribing 
(solid). circle has radius 8 kpc. 
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Fig. 13. The radial dependence of the azimuthally aver- 
aged square of the horizontal magnetic field for the models 
with Ra — 1.5, R^ — 72, a (x lu, but with /,, = 1 (dashed) 
and /,, = 3 (dotted). For comparison the radial profile of 
the azimuthally averaged polarized intensity at A6 cm in 
NGC 1097 is shown with solid line. Note the logarith- 
mic scale on the vertical axes; (B^) = at r = in 
our models. 



Fig. 10. As in Fig. ^, but without any dynamo action: 
Ra = 0, Ruj — 72, /„ = and /,, = 5. In this case the 
magnetic field decays on a time scale of about 10^ yr. 




Fig. 11. Snapshot of typical vectors of the horizontal mag- 
netic field for the positive dynamo number model with 
Ra = —6, Ruj = 36, fa — 0, frf = 5, in a format similar to 
that of Fig. |. 



